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Q\ ', The BRST quantization of the Abelian Proca model is performed using the Batalin- 

Fradkin-Tyutin and the Batalin-Fradkin-Vilkovisky formalism. First, the BFT Hamil- 
Ph! tonian method is applied in order to systematically convert a second class constraint 

system of the model into an effectively first class one by introducing new fields. In find- 
ing the involutive Hamiltonian we adopt a new approach which is more simpler than the 

usual one. We also show that in our model the Dirac brackets of the phase space vari- 
ed ' 

ables in the original second class constraint system are exactly the same as the Poisson 

brackets of the corresponding modified fields in the extended phase space due to the 
linear character of the constraints comparing the Dirac or Faddeev- Jackiw formalisms. 
Then, according to the BFV formalism we obtain that the desired resulting Lagrangian 
preserving BRST symmetry in the standard local gauge fixing procedure naturally in- 
cludes the Stiickelberg scalar related to the explicit gauge symmetry breaking effect 
due to the presence of the mass term. We also analyze the nonstandard nonlocal gauge 
fixing procedure. 
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1 Introduction 

The Dirac method has been widely used in the Hamiltonian formalism 1 to quantize 
the first and the second class constraint systems generally, which do and do not form a 
closed constraint algebra in Poisson brackets, respectively However, since the resulting 
Dirac brackets are generally field-dependent and nonlocal, and have a serious ordering 
problem, the quantization is under unfavorable circumstances because of, essentially, 
the difficulty in finding the canonically conjugate pairs. On the other hand, the quanti- 
zation of first class constraint systems established by Batalin, Fradkin, and Vilkovisky 
(BFV), 2 ' 3 which does not have the previously noted problems of the Dirac method from 
the start, has been well appreciated in a gauge invariant manner with preserving Becci- 
Rouet-Stora-Tyutin (BRST) symmetry. 4 ' 5 After their works, this procedure has been 
generalized to include the second class constraints by Batalin, Fradkin, and Tyutin 
(BFT) 6 ' 7 in the canonical formalism, and applied to various models 8_10 obtaining the 
Wess-Zumino (WZ) actions. n ' 12 

Recently, Banerjee 13 has applied the BFT Hamiltonian method 7 to the second 
class constraint system of the Abelian Chern-Simons (CS) field theory, 14 ~ 16 which 
yields the strongly involutive first class constraint algebra in an extended phase space 
by introducing new fields. As a result, he has obtained a new type of an Abelian WZ 
action, which cannot be obtained in the usual path-integral framework. Very recently, 
we have quantized several interesting models 17 as well as the non-Abelian CS case, 
which yields the weakly involutive first class constraint system originating from the 
non-Abelian nature of the second class constraints of the system, by considering the 
generalized form of the BFT formalism. 18 As shown in all these works, the nature of 
the second class constraint algebra originates from only the symplectic structure of the 
CS term, not due to the local gauge symmetry breaking. Banerjee's and Ghosh 19 have 
also considered the Abelian and non-Abelian (incompletely) Proca model, 20 which 
have the explicit gauge-symmetry breaking term by extending the BFT approach to 
the case of rank-1 non-Abelian conversion compared to the Abelian (rank-0) conversion. 
As a result, the extra field in this approach has identified with the Stiickelberg scalar. 
21 However, all these analysis do not carry out the covariant gauge fixing procedure 



preserving the BRST symmetry based on the BFV formalism hence completing the 
BRST quantization procedure. Furthermore, up to now all above authors do not 
explicitly treat the Dirac brackets in this BFT formalism although the general but 
classical relation between the Dirac bracket and the Poisson bracket in the extended 
phase space are formally reported for the case of Abelian conversion. 7 

In the present paper, the BRST quantization of the Abelian Proca model 20 is 
performed completely by using the usual BFT 7 and the BFV 2 ' 3 formalism. In section 
2, we will apply the usual BFT formalism 7 to the Abelian Proca model in order to 
convert the second class constraint system into a first class one by introducing new 
auxiliary fields. Here, we newly obtain the relation that the well-known Dirac brackets 
between the phase space variables in our starting second class constraint system of 
the Abelian Proca model are the same as the Poisson brackets of the corresponding 
modified ones in the extended phase space without $ — » limiting procedure of the 
general formula of BFT 7 due to essentially the linear character of the constraint. 
It is also compared with the Dirac 1 or Faddeev-Jackiw (FJ) symplectic formalism, 
22 which is to be regarded as the improved version of the Dirac one. Furthermore, 
we adopt a new approach, which is more simpler than the usual one, in finding the 
involutive Hamiltonian by using these modified variables. In section 3, we will consider 
the completion of the BRST quantization based on the BFV formalism. As a result we 
show that by identifying a new auxiliary field with the Stiickelberg scalar we naturally 
derive the Stiickelberg scalar term related to the explicit gauge symmetry breaking 
mass term through a BRST invariant and local gauge fixing procedure according to 
the BFV formalism. We also analyze the nonlocal gauge fixing procedure which has 
been recently studied by several authors. 23 Our conclusions are given in section 4. 



2 BFT Formalism 

2.1 Proca Model and Constraints 

Now, we first apply the usual BFT formalism which assumes the Abelian conversion 
of the second class constraint of the original system 7 to the Abelian Proca model of 



the massive photon in four dimensions, 20 whose dynamics are given by 

S = J d A x [-\f^F^ + \m 2 A»A»\, (1) 

where F^ = d^A u - d u A fM , and g^ = diag(+, -, -, -). 
The canonical momenta of gauge fields are given by 

5S 
tto = -y- « 0, 
SA 

*i = ^~ = F l0 (2) 

SA i v ) 

with the Poisson algebra {A^(x), 7f u (y)} = <5^5(x — y). The weak equality ' ~ ' means 
the equality is not applied before all involved calculations are finished. 1 In contrast, 
the strong equality ' = ' means the equality can be applied at all the steps of the 
calculations. 

Then, Qi = n m is a primary constraint. x The total Hamiltonian is 

H T = H C + f (Pxu^ (3) 



with the multiplier u and the canonical Hamiltonian 

"1 o 1 „ ™, 1 



H r = / d 6 x 



-7T, 



+ -F^F* + -m 2 {(Ay + (A 1 ) 2 } - A n 2 



(4) 



where f2 2 is the Gauss' law constraint, which comes from the time evolution of Qi with 
Ht, defined by 

Q 2 = tfrn + m 2 A° « 0. (5) 

Note that the time evolution of the Gauss' law constraint with Ht generates no more 
additional constraints but only determine the multipler u ~ — diA 1 . As a result, the full 
constraints of this model are Q{ (i,j = 1,2) which satisfy the second class constraint 
algebra as follows 

Aij(x, y) = {fii(x), toj(y)} = -m%£ 3 (x - y), (6) 

detAij(x,y) ^ 0, 

where we denote x = (£, x) and three-space vector x = (x 1 , x 2 , x 3 ) and 612 = —£21 = 1. 



We now introduce new auxiliary fields $* to convert the second class constraints fli 
into first class ones in the extended phase space with the Poisson algebra 

{^(or *„),$*} = 0, 
{&(x),&{y)} = <j j (x,y) = -u ji (y,x). (7) 

Here, the constancy, i.e., the field independence oi uj^(x, y), is considered for simplicity. 
According to the usual BFT method, 7 the modified constraints Qi with the property 

{n i ,f2 i } = 0, (8) 

which is called the Abelian conversion, which is rank-0, of the second class constraint 
(6) are generally given by 

oo 

a(^, tv; &) = a + E $ n) = o; ^ (n) ~ ($T (9) 

fl=l 

satisfying the boundary conditions, fij(A' t ,7r At ;0) = ty. Note that the modified con- 
straints Cli become strongly zero by introducing the auxiliary fields $ 4 , i.e., enlarging 
the phase space, while the original constraints Q{ are weakly zero. As will be shown 
later, essentially due to this property, the result of the Dirac formalism can be easily 
read off from the BFT formalism. The first order correction terms in the infinite series 
7 are simply given by 

n ( j ' ) (x) = Jd 3 yX ij (x,y)&(y), (10) 

and the first class constraint algebra (8) of Cli requires the following relation 

A ij (x,y) + J d 3 w d 3 z X ik (x,w)u kl (w, z)X jl (y, z) = 0. (11) 

However, as was emphasized in Refs. 13, 18, and 19, there is a natural arbitrariness 
in choosing the matrices u/- 7 and X t j from Eqs. (7) and (10), which corresponds to 
canonical transformation in the extended phase space. 6 ' 7 Here we note that Eq. (11) 
can not be considered as the matrix multiplication exactly unless Xji(y, z) is the sym- 
metric matrix, i.e., Xji(y, z) = Xij(z, y) because of the form of the last two product of 
the matrices / d 3 zu kl (uj, z)Xji(y, z) in the right hand side of Eq. (11). Thus, using this 



arbitrariness we can take the simple solutions without any loss of generality, which are 
compatible with Eqs. (7) and (11) as 

^(x,y) = e^ 3 (x-y), 

Xij(x,y) = m(%5 3 (x-y), (12) 

i.e., antisymmetric uj^(x,y) and symmetric Xij(x,y) such that Eq. (11) is the form of 
the matrix multiplication exactly 10 ' 13 ' 17 ^ 19 

Aij(x, y)+ I d 3 ud 3 zX ik (x, u)uo kl (uo, z)X lj (z, y) = 0. 

Note that Xij(x,y) needs not be generally symmetric, while uj l ^{x,y) is always an- 
tisymmetric by definition of Eq. (7). However, the symmetricity of Xij(x,y) is, by 
experience, a powerful property for the solvability of (9) with finite iteration 13 > 17 ~ 19 or 
with infinite regular iterations. 24 

In our model with this proper choice, the modified constraints up to the first order 
iteration term 

Qi = tti + tt^ 

= Qi + m® 1 (13) 

strongly form a first class constraint algebra as follows 

{ik(x) + n?\x), n,-(y) + Qf\y)} = 0. (14) 

Then, the higher order iteration terms 

fiS n+1) = -^*WX*B£> (n > 1) (15) 

n + 2 J 

with 

n n—2 

Bf = £ &f~ m \ Qt\ A ,«) + E {^ n - m) ,fiS m+2) }( $) (16) 

are found to be vanishing without explicit calculation. Here, uik and X k: > are the 
inverse of u lk and X k j, and the Poisson brackets including the subscripts are defined 



by 

dA dB dA dB 



{A,B}( q , p ) 



{A,B} W = £ 



dq dp dp dq ' 

dA dB dA dB 
d& d& d& d& 



where (q,p) and (Q\&) are the conjugate pairs, respectively. 



(17) 



2.2 Physical Variables, First Class Hamiltonian, and Dirac 

Brackets 

Now, corresponding to the original variables A^ 1 and 7fy, the physical variables 
within the Abelian conversion in the extended phase space, A^ and 7r M , which are 
strongly involutive, i.e., 

{n i ,A' i } = 0, {a,^} = 0, (18) 

can be generally found as 

oo 

A' t {A v , n u ; <5> j ) = A" + J2 A^ n \ A^ n) ~ ($ J ') n ; 

n=l 

oo 

*p{A v ,*„;&) = 7r M + 5>(«), 5fW~ (*»")" (19) 

ra=l 

satisfying the boundary conditions, A ,i (A l ',ir l ,',0) = A^ and 7r fi (A u , n u ; 0) = 7r M . Here, 
the first order iteration terms are given by 

A»W = -&u jk X kl {Sl h A»} {A ^ 

m 

TfW = -&uj jk x kl {n l , 7 r lx } {A7r) 

= (m$\0). (20) 

Furthermore, since the modified variables up to the first iterations, A^ + A^ and 
ir v + n^ are found to involutive, i.e., to satisfy Eq. (18), the higher order iteration 



terms 



" n + 

with 



A ^n+l) = —& U J jk X kl (G l *)\ n \ 

T^ 1 ) - —Cf>J,,., Y kl (U ^ n ) 



-VLU jk X kt (HX > (21) 



n n—2 

(cni n) = e i^t m \ a^} m + y: {nt m \^ m+2) } m + R (n+i) , ^ (i) } ($ ), 

m=0 ra=0 

(^)! n) = E {^ n - m) ,^ m) W) + E {^ (n - m) , ^ m+2) } w + R (n+1) , *P}m 

(22) 

are also found to be automatically vanishing. Hence, the physical variables in the 
extended phase space are finally found to be 



A» = 


A" + A^ 1 ' 


7T„ = 


{A° + — $ 2 ,A f + — d^ 1 
m m 


— 


(vro + m^ 1 ,^) 


— 


(fil,7Tj). 



(23) 

Similar to the physical phase space variables A^ and 7f M , all other physical quanti- 
ties, which correspond to the functions of A^ and 7r M , can be also found in principle 
by considering the solutions like as Eq. (19). 7 > 13 > 17 - 19 However, it is expected that 
this procedure of finding the physical quantities may not be simple depending on the 
complexity of the functions. 

In this paper, we consider a new approach using the property 7 ' 25 

K(A»,n^&) = K(A»,n^) (24) 

for the arbitrary function or functional K defined on the original phase space variables 
unless K has the time derivatives: the following relation 

{K(A fi ,ir (1 )M = (25) 



is satisfied for any function K not having the time derivatives because A M and n^ and 
their spatial derivatives already commute with fij at equal times by definition. On the 
other hand, since the solution K of Eq. (24) is unique up to the power of the first class 
constraints fij, 17 ' 19 K(A' 1 ,7f f j i ) can be identified with K^A^, 7r M ; $) modulus the power 
of the first class constraints Qi. However, note that this property is not satisfied when 
the time derivatives exist. 

Using this elegant property we can directly obtain the desired first class Hamiltonian 
Ht corresponding to the total Hamiltonian Ht of Eq. (3) as follows 

Ht(A»,7t u ;<5> 1 ) = H T (A^,n v ) 

= H T (A", n v ) + [ d 3 x Wm 1 ) 2 + i($ 2 ) 2 + -dA^^i - -$ 2 ^ 2 • 
J \_2 2 m m 

(26) 

On the other hand, since the first class Hamiltonian H c corresponding to the canonical 
Hamiltonian H c of Eq. (4) can be similarly obtained as follows 



H c (A»,ir v ;&) = H c (A»,n u ) 

= H c (A»,ir u ) + Jd 3 x 



-(di® 1 ) 2 + i($ 2 ) 2 + m&&A i - -$ 2 fi 2 

2 2 m 



(27) 



Eq. (26) can be re-expressed as follows 

Et{A\ 7T„; $<) = H C (A", 7r„; $*) - J d^A^ (28) 

as it should be according to Eq. (3). 

This is the same result as the usual approach (See the Appendix A) and, by con- 
struction, both Ht and H c are automatically strongly involutive, 

{n h Hr} = 0, (29) 

{Q t ,H c } = 0. 

Note that all our constraints have already this property, i.e., f2j(A M , 7r M ; $) = f2j(A M , 7f M ). 
In this way, the second class constraints system fij(A^, 71^) ps is converted into the 
first class constraints one fi^A^, 7T M ; $) = with the boundary conditions fij|$ = o = O^. 



On the other hand, in the Dirac formalism 1 one can make the second class con- 
straint system Qj rs into the first class constraint one Qi(A^, 7r M ) = only by deform- 
ing the phase space (A^, n^) without introducing any new fields. Hence, it seems that 
these two formalisms are drastically different ones. However, remarkably the Dirac 
formalism can be easily read off from the usual BFT-formalism 7 by noting that the 
Poisson bracket in the extended phase space with $ — ► limit becomes 

{A,B}\„=o = {A,B}-{A,n k }A kk '{n k ,,B} 

= {A,B} D (30) 

where A kk = —X lk ujwX lk is the inverse of A^/ in Eq. (6). About this remarkable 
relation, we note that this is essentially due to the Abelian conversion method of 
the original second class constraint: In this case the Poisson brackets between the 
constraints and the other things in the extended phase space are already strongly zero 

{n t ,A} = o, 

{&,&,} = 0, (31) 

which resembles the property of the Dirac bracket in the non-extended phase space 

{n h A} D = o, 

{n t ,n 3 }D = o, (32) 

such that 

{f2i,A}|*=o = {n l ,A}* = o, 

{^,^}| $=0 = {n i ,n j y = o (33) 

are satisfied for some bracket in the non-extended phase space { , }*. However, due 
to the uniqueness of the Dirac bracket 26 it is natural to expect the previous result (29) 
is satisfied, i.e., 

{ , }• = { , }d (34) 



10 



without explicit manipulation. Moreover we add that, due to similar reason, some 
non-Abelian generalization of the Abelian conversion as 

{n t ,A} = aij & + a ijk &$ k + ---, 
{f^fy} = ^ fc $ fc + ^$ fe $' + ... (35) 

also gives the same result (30) with the functions «„, otijk, /%&, etc, • • • of the original 
phase variables A^, ix v . As an specific example, let us consider the brackets between 
the phase space variables in Eq. (23). The results are as follows 

{A (x),Ai(y)}\* =0 = {A°(x),Ai(y)} = ±d> x 6(x-y), 

{A (x),A°(y)}\* =0 = {A°(x),A°(y)} = 0, 

{A^x),A k (y)}\ <s>=0 = {A\x),A k {y)} = 0, 

{n M (x),n,( y )}| $=0 = {n„(a;),n„(j/)} = 0, 

{A^ILjiy)}^ = {A\x)^{y)} = M(* " v), 

{A (x),n v (y)}|*=o = {A°(x),UM} = 0, 

{#(x),no(s/)}|*= = K(x),n (y)} = (36) 

due to the linear correction (i.e., only the first order correction) of the modified fields 
A^ and 7r M , which are the same as the usual Dirac brackets. 19 Hence, in the case of 
the phase space variables of the model the well-known Dirac brackets of the fields are 
exactly the Poisson brackets of the corresponding modified fields in the extended phase 
space contrast to the general formula (30). Note that the FJ symplectic formalism, 
22 which is the improved version of the Dirac method, also gives the same result (See 
Appendix B for this matter). 

Now, since in the Hamiltonian formalism the first class constraint system without 
the CS term 13 ' 18 indicates the presence of a local symmetry, this completes the op- 
eratorial conversion of the original second class system with the Hamiltonian H? and 
the constraints fii into first class one with the Hamiltonian Ht and the constraints fii. 
From Eqs. (13) and (28), one can easily see that the original second class constraint 
system is converted into the effectively first class one if one introduces two fields, which 
are conjugated with each other in the extended phase space. Note that for the Proca 

11 



case the origin of the second class constraint is due to the explicit gauge symmetry 
breaking term in the action (1). 



2.3 Corresponding First Class Lagrangian 

Next, we consider the partition function of the model in order to present the La- 
grangian corresponding to Ht in the canonical Hamiltonian formalism. However, the 
result is the same with H c . As a result, we will unravel the correspondence of the 
Hamiltonian approach with the well-known Stiickelberg's formalism. First, let us iden- 
tify the new variables $* as a canonically conjugate pair (p, ix p ) in the Hamiltonian 
formalism, 

(&)^( mp ,±-n p ), (37) 

satisfying Eqs. (7) and (12). Then, the starting phase space partition function is given 
by the Faddeev formula 3 ' 27 as follows 

r 2 

Z= VA^V^VpVnp [] S(Qi)6(T j)det | {J^r,} | e lS , (38) 

where 



S = J (Tx [n p A^ + n p p - ri T ) , (39) 

with the Hamiltonian density 7Yt corresponding to the Hamiltonian Ht of Eq. (25), 
which is now expressed in terms of (p, n p ) instead of $\ Note that the gauge fixing 
conditions Tj are chosen so that the determinant occurring in the functional measure 
is nonvanishing. Moreover, Tj may be assumed to be independent of the momenta so 
that these are considered as the Faddeev-Popov type gauge conditions. 10 . 13 . 17 . 27 . 28 

Before performing the momentum integrations to obtain the partition function in 
the configuration space, it seems appropriate to comment on the involutive Hamilto- 
nian. If we directly use the above Hamiltonian following the previous analysis done 
by Banerjee et ai, 19 we will finally obtain the non-local action corresponding to this 
Hamiltonian due to the existence of (9Vj) 2 -term in the action when we carry out the 
functional integration over n p later. Furthermore, if we use the above Hamiltonian, we 
can not also naturally generate the first class Gauss' law constraint Q2 from the time 

12 



evolution of the primary constraint fii, i.e., {Cl\, H T } = 0. Therefore, in order to avoid 
these unwanted situations, we use the equivalent first class Hamiltonian without any 
loss of generality, which differs from the involutive Hamiltonian (26) by adding a term 
proportional to the first class constraint Q 2 as follows 

Wt =Ht+ I d 3 x^Q 2 . (40) 

J m z 

Then, we have the desired first constraint system such that 

{n^Hr'} = n 2 , 

{n 2 ,Hr'} = 0. (4i) 

Note that when we operate this modified Hamiltonian on physical states, the difference 
is trivial because such states are annihilated by the first class constraints. Similarly, 
the equations of motion for gauge invariant variables will also be unaffected by this 
difference since Q 2 can be regarded as the generator of the gauge transformations. 
Now, we consider the following effective phase space partition function 

r 2 

Z = / Vit li VA"V'n p Vp J] S(Q i )5(T j )det | {Sl^V,} | e lS \ 

S' = f d 4 x (tt m > + n p p - H' T ). (42) 

The 7To integral is performed trivially by exploiting the delta function 8{£l\) = 5(tyq + 
m 2 p) in Eq. (42). On the other hand, the other delta function 5(Q 2 ) = 5(d l TT i + m 2 A° + 
Up) can be expressed by its Fourier transform with Fourier variable £ as follows 

6{n 2 )= [V£e- i f d4x &z. (43) 

Making a change of variable A — > A + £, we obtain the action 

S = J d A x [n t A* - m 2 p(A* + £) + n p p - ^ - ±i^F*" + ^m 2 (A ) 2 + \m 2 A % X 

+ AWm - m 2 d^p " -j^ + \m*d iP &p - ^ p - \m 2 e\ , (44) 

where the corresponding measure is given by 

[Dp] = V^Vn^A^VpVTr.YliSlTjiA + £, A\ tt 4 , p)]}det | {Sl u T,} | . (45) 
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Performing the Gaussian integral over 7Tj, this yields the intermediate action as follows 

S u = Jd 4 x[~F^F^ + ±m 2 A p A» 

+ n p (p - e - 2^tt p ) - m 2 p(i° + - m 2 d % A*p + \m 2 d lP &p - ^m 2 £ 2 ].(46) 

To realize the Stiickelberg term through the BFT analysis, we choose the Faddeev- 
Popov-like gauges 10 ' 13 ' 17 ' 27 ' 28 j which do not involve the momenta. After the Gaussian 
integration over 7r p , we finally obtain the well-known action up to the total divergence 
by identifying the extra field p with the Stiickelberg scalar as follows 

S = J d 4 x[-^F, u F^ + \m 2 {A, + d,p) 2 ], (47) 

which is invariant under the gauge transformations as 5A P = <9 M A and dp = —A. As 
expected, the Stiickelberg scalar p is introduced in the mass term. 

It seems to appropriate to comment on the original unitary gauge fixing. If we 
choose the gauge as follows 

Ti = (p, *„), (48) 

we get 

{^(x),r j (y)} = e,^ 3 (x-y). (49) 

Then, integrating over the variables p and ir p we reproduce the original partition func- 
tion as follows 

Z = fvA^S d ' x [-aW'-X^"]. (50) 

The physical meaning of this result is that the original action can be regarded as a 
gauge-fixed version of the first constraint system (18) and (40). Note that this gauge 
fixing is consistent because when we take the gauge fixing condition p ~ 0, the condition 
T\ p ~ is naturally generated from the time evolution of p, i.e., p = {p, H u } = — ^n p m 
0, where the Hamiltonian H u corresponds to the intermediate action S u . 



3 BFV-BRST Gauge Fixing 
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3.1 Basic Structure 

In this subsection, we briefly recapitulate the BFV formalism 2 ' 3 which is applicable 
for general theories with the first-class constraints. For simplicity, this formalism is 
restricted to a finite number of phase space variables. This makes the discussion simpler 
and conclusions more apparent. 

First of all, consider a phase space of the bosonic canonical variables q l , Pi (i = 1, 2, 
• • •, n) in terms of which canonical Hamiltonian H c (q\pi) and constraints Q a (q l ,Pi) ~ 
(a = 1, 2, • • •, m), which being bosonic also, are given. We assume that the constraints 
satisfy the following constraint algebra 2 ' 3 ' 9 

{Q a , Q b } = U^ b fi c , 

{H c ,n a } = v£n b , (5i) 

where the structure coefficients U^ b and V£ are functions of the canonical variables. 
We also assume that the constraints are irreducible, which means that there locally 
exists an invertible change of the variables such that Q a can be identified with the 
m-unphysical momenta. 

In order to single out the physical variables, we can introduce the additional bosonic 
conditions T a (q l ,pi) ps with det |{T a , Q b }\ ^ at least in the vicinity of the constraint 
surface r a ps and Q a rs 0. Then, T a play the roles of gauge- fixing functions. That is 
to say, from the condition of time stability of the constraints, a family of phase space 
trajectories is possible. By selecting one of these trajectories through the conditions 
of T a ps 0, we can get the 2(n — m) dimensional physical sub-phase space denoted by 
q*, p*. And then, r a (q l ,pi) can be identified with the m-unphysical coordinates. 

The described dynamical system with the partition function 

z = J[dq l d Pt } 5(n a )s(r b ) det |{r 6 , n a }\ e ^-w-^) (52) 

is then completely equivalent to the effective quantum theory with the following par- 
tition function 

Z eq . = J[dq*dp*] e '/*W-iW(«V)] (53) 

which only depending on the canonical variables q*,p* of the physical sub-phase space, 

15 



And the constraints fi a ~ and r a w together with the Hamilton equations may 
be obtained from a action 



s = Jdt ( Pt q i -h c - x a n a + 7r a r a ) 



(54) 



where A a and ir a are the bosonic Lagrange multiplier fields canonically conjugated to 
each other, obeying the Poisson bracket relations 



{A a , 7T a } = 81 . 



Note that the gauge-fixing conditions contain A a in the following general form 



(55) 



T a = \ a + x a (q\ Pl ,X a ), 



(56) 



where x a are arbitrary functions. And we can see that the Lagrange multiplier A a be- 
come dynamically active, and n a serve as their conjugate momenta. This consideration 
naturally leads to the canonical formalism in an extended phase space. 

In order to make the equivalence to the initial theory with the constraints in the 
reduced phase space, we may introduce two sets of canonically conjugate, fermionic 
ghost coordinates and momenta C a , V a and V a , C a such that 



{C a ,V b } = {V a ,Cb} = s a b 



(57) 



with the super-Poisson bracket 



{A,B} 



SA 



5_B_ 

Sp 



-(-1) 



VA'IB 



8_B_ 



SA 

r S P 



where tja denotes the number of fermions called ghost number in A, and subscript "r" 
and "/" right and left derivatives. 

The quantum theory is now defined by the extended phase space functional integral 



Z* = [dq l d Pi }[d\ a dn a }ldC a dV a }[dV a dC a }e lb * 



(58) 



where the action is 



S* = / dt{ Pl q l + n a X a + V a C a + C a V a -H m + {Q, ^}}. 



(59) 
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Here, the BRST-charge Q and the fermionic gauge-fixing function \1/ are defined by 

Q = c a n a - l -c b c c u%v a + v a n a , 

tf = C aX a + V a \ a , (60) 

respectively. H m is the BRST invariant Hamiltonian, called the minimal Hamiltonian, 

H m = H c + C a V h a V b . (61) 

The measure in Zq is the Liouville measure on the covariant phase space. Further- 
more, if we choose the fermionic gauge-fixing function \1/ properly, 3 ' n we can obtain 
manifestly covariant expression. And the equivalence of the dimensionality 2n + 6m 
in the extended phase space, including the canonical ghost variables, to the original 
dimensionality In — 2m in the reduced phase space can be seen by identifying the ghost 
variables with the negative-dimensional canonical degree of freedom, which is suggested 
by the Parisi-Sourlas' original work related with the superrotation Osp(l,l|2) in the 
extended phase space. 29 

In order to derive the BRST gauge-fixed covariant action for the Abelian Proca 
theory, according to the above BFV formalism in the extended phase space, let us 
introduce the ghosts and anti-ghosts together with auxiliary fields as follows 

(C*,P t ), (V'A), (N\B t ), (62) 

where % — 1, 2. The nilpotent BRST-charge Q, the fermionic gauge-fixing function \I/ 
and the minimal Hamiltonian H m in our case are 

Q = fdx [C% + T^Bi], 

* = Jdx[C lX l + V t N l ], 

H m = Ht' - jdx [V 2 C l ], (63) 

where x 1 — ^4°> X 2 = diA 1 + ^-82, and a is an arbitrary parameter. 

The BRST-charge Q, the fermionic gauge-fixing function \t>, and the minimal Hamil- 
tonian H m satisfy the following relations, 

{Q,H m } = 0, 

Q 2 = {Q,Q} = 0, 
{{^,Q},Q} = 0, (64) 
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which being the conditions of physical subspace after the operator quantization 

'A,B\ =ih{A,B} (65) 

for the quantum operators A and B corresponding to the classical functions A and B 
when there is no operator ordering problem. 
The effective action is 

S eff = Jd 4 x [ 7T A° + itiA 1 + n p p + B 2 N 2 + V % C + C 2 V 2 ] - H total (66) 

where H total = H m + {Q,^}, and also / d^x^N 1 + C 1 V 1 ) = {Q,Jd 4 x CN 1 } terms 
are suppressed by replacing \ l with x 1 + N 1 just like the cases in Refs. 8, 9. 



3.2 Local Effective Action 

Now, in order to derive the covariant effective action we first perform the path 
integration over the fields Bi, N 1 , C\,V l , V\, C 1 , A , and tcq by using of the Gaussian 
integration. Then, we obtain 

S eff = fd 4 x [ TiiA 1 + tt p p + BN + VC + CV 

- -(vr,) 2 - -^(n p ) 2 - -FijF** - -m 2 (A 1 ) 2 
2 K J 2m 2 v pJ 4 3 2 K J 

- m 2 pd % A % + ]^m 2 d iP d l p - N{d l n t + n p ) 

- B{d t A l + -aB) ~ d£&C - VV } (67) 

with N 2 = N, B 2 = B, C 2 = C, C 2 = C, V 2 = V, and V 2 = V. Using the variations 
over re 1 , 7T p , V and V, we obtain the following relations 



m = A 1 - d t N, 

7r p = m 2 (p-N), 

V = -t, V = C, (68) 

and identifying with N = —A , we get the usual local form of the covariant effective 
action 

S eff = jd 4 x [-\f^F» v + l -m 2 (A p + 3,p) 2 + A^B - ^a(B) 2 - dfid'C ], (69) 



which is invariant under the standard BRST transformation 

5 B A„ = -Xd^C, S B p = XC, 
S B C = 0, 5 B C = -XB, 5 B B = 0, (70) 

which is local and covariant one. This completes the procedure of BRST invariant 
(here standard) gauge fixing with the local action according to the BFV formalism. 

Therefore, we see that the auxiliary BF field p is exactly the well-known Stiickelberg 
scalar in Eq. (69). 



3.3 Nonlocal Effective Action 

Although in the previous subsection we have performed the integration over V and 
V but not over C and C, it is not impossible to consider the opposite procedure, i.e., 
the integration over C and C but not over V and V which is dual to the previous 
integration. To this end, we consider the BFV formalism in the previous subsection 
up to the point where the integration over the momentum 7r p was performed and the 
following effective action was obtained: 



S eff = jd A x[-- A F, u F^+ l -m\A, + d,pf 



-A^d^B + \a{Bf - diCd l C + VC + CV-VV]. (71) 

This action is invariant under the BRST transformation which have the form 

5 B A = -XV, S B A t = -XdiC, 5 B p = XC, 
5 B C = 0, 5 B C = -XB, 5 B B = 0, 
S B V = 0, 5 B V = -X(-d t F 0l + m 2 (p + A )). (72) 

Now, let us perform the integration over C, C instead of their conjugated ones V, 
V. First, performing the integration over the ghost field C, we get the following delta 
function 

Sfr&C-V) = fet{d i &)8{C- 1 j^P). (73) 

19 



Next, performing the integration over C, we get the unusual non-local form of the 
effective action 

Self = jd'x i-±F, u F^ + ±m\A, + d,pf-A^B + ±a(B) 2 + V^V-VV] (74) 

which is non-covariant. Notice that the appearance of the nonlocal term in the ghost 
action was a result of the unusual integration. But, we may find that this form is also 
obtained by the change of variables 

c - w v - z -* T (75) 

in the Eq. (69). Under these replacements, we have the nonlocal BRST charge Q' 
given by 

Q' = Jd 3 x[BC + (-d t F * + m 2 (p + A°))J^V}. (76) 

Then, the effective action is invariant under the following nonstandard BRST trans- 
formation as 

S' B V = 0, 5' B V = -XB, 5' B B = 0, (77) 

which is non-local and non-covariant and hence can be categorized as the symmetry 
recently proposed in QED by Lavelle and McMullan. 23 Moreover, this nonlocal BRST 
symmetry yields a conserved current through the Noether's theorem as follows 

f» = W^* + ™'(A. + M^ + B6,±V. (78) 

This completes the procedure of the BRST invariant (here nonstandard) gauge fixing 
with the non-local action according to BFV formalism. As a result, we have recognized 
that these nonlocal symmetry and the conserved current are nothing but the original 
local theory performing the change of variable (75). 



4 Conclusion 

In conclusion, we have applied the BFT and the BFV method to covariantly quan- 
tize the second class constraint system of the Abelian Proca model without spoiling 
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the unitarity. First, by applying the BFT method, we have systematically converted 
the second class system of the model into the effectively first class one in the extended 
phase space. We have shown the relation that, due to the linear character of the con- 
straint, the Dirac brackets between the phase space variables in the original second 
class system are exactly the Poisson brackets of the corresponding modified ones in the 
extended phase space without $ — > limiting procedure of the general formula of BFT 
7 by comparing the Dirac or Faddeev-Jackiw symplectic formalism. 

Furthermore, we have noted that, like as this relation, the general relation (30) of 
the Dirac brackets in the non-extended phase space and the Poisson brackets in the 
extended phase space is essentially due to the Abelian conversion (31) of the second 
class constraint into the first class one and we have added also that some more general 
conversion method (35) may be considered without spoiling this nice relation. 

Moreover, we have adopted a new approach, which is more simpler than the usual 
one, in finding the involutive Hamiltonian by using these modified variables according 
to the important property (24). Now, with this first class constraint we have applied 
the BFV method to covariantly quantize the Proca model without spoiling the unitar- 
ity. By identifying a new auxiliary field with the Stuckelberg scalar we have naturally 
derived the Stuckelberg scalar term related to the explicit gauge symmetry breaking 
mass term through the usual local gauge fixing procedure with the standard BRST 
symmetry according to the BFV formalism. We have also analyzed the nonlocal gauge 
fixing procedure with the nonstandard BRST symmetry which has been recently stud- 
ied by several authors. 23 

As final remarks, we first note that although we have successfully applied the rather 
simple Abelian case, it is not clear whether the non- Abelian generalization of our model 
is possible or not in a priori. However, considering the recent failure of the complete 
conversion of the second class constraints of this non- Abelian model into the first class 
ones, 19 which does not use the symmetric Xij(x, y)-matrix, and the power of the 
formalism with the symmetric Xij(x, y) as noted in Section 2.1 we expect that whether 
we can find the symmetric Xij(x, y) or not is crucial point for the successful application 
of our BFT formalism. We are in progress in this direction. In another direction, a 
new formalism may be considered to solve this problem, but it is not clear whether the 
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complete conversion of the system into the first class one is possible and furthermore 
that formalism is equivalent to the BFT formalism or not. 30 On the other hand, it 
is interesting to note that the non-Abelian model of the Chern-Simons theory allows 
the symmetric Xij(x,y) and the solution of the first class system is found by finite 
iterations. 18 



Appendix A 

In this appendix, we derive the first class Hamiltonians (26) and (27) in the ex- 
tended phase space corresponding to the total Hamiltonian Ht of (3) and canonical 
Hamiltonian H c of (4) by using the usual straightforward approach. 13 > 17 - 19 > 28 Let us 
first consider the first class Hamiltonian Ht corresponding to Ht- It is given by the 
infinite series, 

oo 

H T = H T + Y J H^ ) ] 4 n) ~(c&r, (79) 

n=l 

satisfying the initial condition, Ht^^, A^; 0) = Ht- The general solution 7 for the 
involution of Ht is given by 

H { t ] = — I d 3 xd 3 yd 3 z ^(x^jix, y)X ]k {y, z)G { £~ l \z) (n > 1), (80) 

where the generating functions G£ are given by 

Gf = {<l?\H T }, 

Gf = {Sll \HF>} o + {Ql 1 \H!r 1) }o (n>l), (81) 

where the symbol O in Eq. (80) represents that the Poisson brackets are calculated 
among the original variables, i.e., O = (A M ,7r M ). Here, LVy and X u are the inverse 
matrices of a/ 7 and X^ respectively as in the text. Explicit calculations yields 

Gf = Sh, 

Gf = m^i, (82) 

which are substituted in (80) to obtain H T , 

H T l) = I d 3 x l—tdidiP 1 )^ - — $ 2 fi 2 l . (83) 

J Ym m 

22 



This is inserted back in Eq. (81) in order to deduce G\ as follows 



G? 
GP 



m$ 2 , 
mdidi^ 1 . 



r(2) 



(i) 



Then, we obtain Hj, by substituting G\ in Eq. (76) 



H. 



(2) 



drx 



-(ft* 1 )^* 1 ) - \{® 2 ) 2 



J4) 



55) 



6) 



Finally, since 

Gf } = (n>2), 

due to the proper choice (12) we obtain the complete form of the Hamiltonian H as 
follows 

H t = H t + h£ ) +H¥\ (87) 

which, by construction, is strongly involutive, 



{Qi,H T } = 0. (88) 

Similarly, for the canonical Hamiltonian we can easily obtain it's first class Hamil- 
tonian H r as follows 



h c = h c + h^ + hP, 



where 



tfW 



d x 



m&idiA* 
1 



1 
m 



$ 2 a 



H^ = J d\[^){d^) - '-{^f 



(89) 

(90) 
(91) 



Here, we used 



Gf 
G?> 

G\ (n) 



ft 2 , 

m 2 <9^, 

m<l> 2 , 

— mdidi^ 1 . 

(ra>2). 



(92) 
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Note the differences in G 2 , G 2 and ^-dependent terms in H^ 1 ' and H^ 2 > for the total 
and canonical Hamiltonians. Moreover, H c is also, by construction, strongly involutive, 

{Qi, H c } = 0. (93) 



Appendix B 

In this appendix, we obtain the FJ symplectic brackets comparing them with both 
the orthodox Dirac brackets and the modified Poisson brackets (36) in the extended 
phase space in the section 2.2. 

According to the FJ formalism, 22 which is to be regarded as the improved version 
of the Dirac one, we rewrite the first order Lagrangian corresponding to the Proca 
model (1) as 

£(°) = n t A l -H {0) , (94) 

where the conjugate momenta (2) of the gauge fields and the canonical Hamiltonian 
7i c in which we denote it as 7i^ showing the iterative nature of the formalism are 
used. 

In order to find the FJ symplectic brackets we introduce the sets of the symplectic 
variables ^°' k and the conjugate momenta a k as follows 

4 0) = (0,7r„0), (95) 

which are usually read off from the form of the canonical sector of the first order 

Lagrangian (94), respectively. 

Then, the dynamics of the model is governed by the invertible symplectic two-form 

matrix such that 

f (., (jj) *■!%) *,%) 

through the equations of motion 

fM = /rfV<n^)f^g, 07) 
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where f^°' kl (x, y) is an inverse of /^ (x, y). However, in the Proca model the symplectic 
two-form matrix is given by 

/ \ 

fif(x,y)=\ -Sy \S 3 (x-y) (98) 

\o Sa o J 

showing the matrix /^ (x, y) is singular. As it happens, the symplectic two-form matrix 
has a zero mode, i.e., v\ = (i>i,0,0), where V\ is an arbitrary function. Furthermore, 
this zero mode generates constraint f^ 1 ) in the context of the FJ formalism 22 as follows 

= J d 3 x vi(x)jj— -J 'd 3 yH(°\y) 

d 3 xv 1 {d l n t + m 2 A°) 

d 3 x v^, (99) 

which will be added into the canonical sector of the Lagrangian (94) enlarging the 
symplectic phase space with the Lagrange multiplier a. Then, the iterated, first order 
Lagrangian is given by 

£(!) = ni A l + (8% + m 2 A°)a - H (1 \ (100) 

where the corresponding first iterated Hamiltonian liS 1 ^ is given by 

W (1) (0 lod)=o = \^f + \F^ + \m\A'Y + \m\A*y. (101) 

The situations we stand is exactly the same as before except we now have the first 
order Lagrangian (100) and the Hamiltonian (101). In other words, we can set again 
the symplectic variables and the conjugate momenta as follows 

e )k = (A°,A\n t ,a), 
4 1} = (O^CdVi + mM ) (102) 

reading off from the Lagrangian (100). From this set of the variables, the first iterated 
symplectic two-form matrix is given by 



/£W) 



/ 











-Sij 


m 2 \ 



V 



— m 


Sij 





d 1 


/ 



S 3 (x-y), (103) 
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and its inverse matrix is easily obtained 

/ -\dt -4j \ 

m A x m* 



r w (x,y) 



-^di o St, o 
-Sij 

V A ooo/ 



S 3 (x-y). (104) 



Now, according to the FJ formalism, this inverse symplectic two-form matrix gives the 
symplective brackets of the Proca model 

{e(x),e(y)}sym P = f kl{1) (x,y) (105) 

in the case of having the invertible symplectic matrix, i.e., at the final stage of iteration, 
such that 



[X), A {yjfsymp 


= —j d i s ( x -y) 


yx), A yy) jsymp 


= 0, 


\X), A [y) Jsymp 


= 0, 


fj,\X) : i-i-uyy ) J symp 


= 0, 


\X)i 11 j{y ) Jsymp - 


= 8ij8(x-y), 


yx ) , ii;, yy ) j symp 


= 0, 


(x), n (?/) Jsymp z 


= 



(106) 

showing that the symplectic brackets are exactly same both as the Dirac brackets and 
the modified Poisson brackets in the extended phase space (36). 
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